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^ ■ ABSTRACT 

' We derive the general resummed effective Lagrangian for Higgs-mediated flavour-clianging neutral- 

5— ( ' current (FCNC) interactions in the Minimal Super symmetric Standard Model (MSSM), without 



resorting to particular assumptions that rely on the squark-mass or the quark- Yukawa structure 
of the theory. In our derivation we also include the possibility of explicit CP violation through 
the Cabibbo-Kobayashi-Maskawa mixing matrix and soft supersymmetry-breaking mass terms. 
The advantages of our resummed FCNC effective Lagrangian are explicitly demonstrated within 
the context of phenomenologically motivated scenarios. We obtain new testable predictions in 
the large tan /? regime of the theory for CP-conserving and CP-violating observables related to 
the K- and i?-meson systems, such as AMk,b, ^k-, e'/^) ^{Bd,s and their associated 

leptonic CP asymmetries. Finally, based on our resummed FCNC effective Lagrangian, we can 
identify configurations in the soft supersymmetry-breaking parameter space, for which a kind 
of a Glashow-Iliopoulos-Maiani-cancellation mechanism becomes operative and hence all Higgs- 
mediated, tan /^-enhanced effects on K- and i?-meson FCNC observables vanish. 
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1 Introduction 



The appearance of too large flavour-changing neutral-current (FCNC) interactions of Higgs 
bosons to fermionic matter is a generic feature of SU(2)LxU(l)y theories with two and 
more Higgs doublets. Unless there is a symmetry to forbid these Higgs-mediated FCNC 
interactions to occur in the bare Lagrangian of the model their unsuppressed existence 
will inevitably lead to predictions for rare processes in the kaon and 5-meson systems 
that violate experimental limits by several orders of magnitude § . In the minimal 
realization of softly-broken supersymmetry (SUSY), the Minimal Supersymmetric Standard 
Model (MSSM), the holomorphicity of the superpotential prevents the occurrence of Higgs- 
boson FCNCs by coupling the one Higgs-doublet superfield, Hi, to the down-quark sector, 
and the other one, H2, to the up-quark sector. However, the above holomorphic property of 
the superpotential is violated by finite radiative (threshold) corrections due to soft SUSY- 
breaking interactions 0, |[. As a consequence, Higgs-mediated FCNCs reappear at the 
one-loop level, but are naturally suppressed for low and intermediate values of tan/? = 
{H2)/{Hi), i.e. for tan/3 < 20. For larger values of tan/5, e.g. tan/5 > 30, the FCNCs 
partially overcome the loop suppression factor 1/(167?^) and become phenomenologically 
relevant especially for the K- and 5-meson systems. 

Recently, the topic of Higgs-boson FCNCs in the large-tan /5 limit of the MSSM has 
received much attention [|,|^,|^,||,miy,|Tl|,|12|,|l^,[l^,|15|,0]. Several approaches have been 
devised to implement the non-holomorphic finite radiative corrections into the phenomeno- 
logical analysis of FCNC processes, such as K^K^ and B^B^ mixings, B and 
Bs — ^ In most cases, however, the suggested approaches to threshold radiative ef- 

fects involve certain explicit or implicit assumptions pertinent to the squark-mass and the 
quark- Yukawa structures of the theory, such as the dominance of the top quark in the FCNC 
transition amplitudes. We term the latter assumption the t-quark dominance hypothesis. 
On the other hand, some of the approaches neglect higher-order terms in the resumma- 
tion of threshold corrections to (i-quark Yukawa couplings, which become important in the 
large-tan (3 regime of the theory. 

In this paper we derive the effective Lagrangian that properly takes into account the 
resummation of higher-order threshold effects on Higgs-boson FCNC interactions to down 
quarks. To accomplish this in Section 2, we avoid the imposition of particular assump- 
tions on the structure of the soft squark masses and the quark- Yukawa couplings of the 
theory. Moreover, we do not rely on specific kinematic approximations to the transition 
amplitudes, such as the aforementioned t-quark dominance hypothesis in the FCNC matrix 
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elements. In our derivation of the effective Lagrangian, we also consider the possibility of 
CP violation through two sources: (i) the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix [|I^ and (ii) the soft SUSY-violating mass terms. As we explicitly demonstrate 
in Section 3, our resummed FCNC effective Lagrangian gives rise to new testable predic- 
tions for CP-conserving as well as CP-violating observables related to the K- and 5-meson 
systems. In the same section, we qualitatively discuss the implications of tan /5-enhanced 
Higgs-mediated interactions for the direct CP- violation parameter e'/ e in the kaon system. 
Section 4 is devoted to our numerical analysis of a number of K- and 5-meson observables, 
such as AMk, €k, ^Mb^^, B{Bd,s and their associated leptonic CP asymme- 

tries. In particular, based on our resummed FCNC effective Lagrangian, we are able to 
identify configurations in the soft SUSY-breaking parameter space, for which a kind of a 
Glashow-Iliopoulos-Maiani-cancellation mechanism (GIM) |18[ becomes operative in the 
Higgs-ci-quark sector. As a result, all Higgs-mediated, tan /3-enhanced effects on K- and 
5-meson FCNC observables vanish. Finally, Section 5 summarizes our conclusions. 



2 Resummed FCNC effective Lagrangian 

In this section, we derive the general form for the effective Lagrangian of Higgs-mediated 
FCNC interactions in the CP-violating MSSM. For this purpose, we also consistently resum 
the tan /3-enhanced radiative effects on the cZ-quark Yukawa couplings 0. First, we analyze 
a simple soft SUSY-breaking model based on the assumption of minimal flavour violation [^, 
|lO|,|T3|,0], where the CKM matrix is the only source of flavour and CP violations. We find 
that even within this minimal framework, the usually neglected c-quark contribution to 
Higgs-mediated FCNC interactions may be competitive to the t-quark one in certain regions 
of the parameter space. After having gained some insight from the above considerations, 
we then extend our resummed effective Lagrangian approach to more general cases that 
include a non-universal or hierarchical squark sector as well as CP violation originating 
from the CKM matrix and the soft SUSY-breaking parameters. 

Before discussing the most general case, let us first consider the following simple form 
for the effective Yukawa Lagrangian governing the Higgs-mediated FCNC interactions in 
the quark sector 



— Cy — d^nii-d 



< + (Eg + E„hl,h„)] 4 + $°4h,< + H.c, (2.1) 



where $5 2 electrically neutral dynamical degrees of freedom of the two Higgs 
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Figure 1: Non-holomorphic radiative effects on the d-quark Yukawa couplings induced by 
(a) gluinos g and (b) charged Higgsinos hi2- 



doubletsQand the superscript '0' on the d- and w-type quarks denotes fields in the interaction 
basis. In ( p.l| ), hd and hu are 3 x 3-dimensional down- and up-quark Yukawa matrices, 
and SBi 



E, 



1 



**r/2 2 I |2\ 



167r2 



(2.2) 
(2.3) 



are finite non-holomorphic radiative effects induced by the diagrams shown in Fig. |l]. In 
the above, the loop integral I{x, y, z) is given by 

xylnixly) + yzhiiijlz) + xz\\i{z/x) 



I{x,y,z) 



(2.4) 



{x -y){y- z) {x- z) 

with /(x, x,x) = l/(2x). To keep things simple in the beginning, we assume that in ( |2.2| ) 



and 



the bilinear soft SUSY-breaking masses of the squarks, m. 



and the 



trilinear soft Yukawa couplings Au = A^^d are flavour-diagonal and universal at the soft 
SUSY-breaking scale Msusy- We also neglect the left-right mixing terms ul-ur and di-dR 
in the squark mass matrices. The consequences of relaxing the above assumptions will be 
discussed later on. 

From (|2.1|) , we can easily write down the effective Lagrangian relevant to the effective 
d- and u-type quark masses: 



1 + tan /3 E 



Euh|,h, 



dl 



^4h,< + H.c. (2.5) 



Our next step is to redefine the quark fields as follows: 



U^UL. dl = U^VdL, 



Ui 



U^Ur, 



4 



dR 



[2.6) 



^Throughout the paper, we follow the notation and the CP-phase conventions of 
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where U^, U^, and V are 3-by-3 unitary matrices that relate the weak to mass eigen- 
states of quarks. Evidently, V is by construction the physical CKM matrix. Substitut- 
ing ( pl6| ) into ( pISD yields 



Vl 



~ -^mass — dR hjW^ 



1 + tan/5 (^Eg + Eu|hu| 



V2 



Y di + —^URhuUL + H.c. 
v2 

(2.7) 



where Md and Mu are the physical d- and w-quark mass matrices, respectively. Consistency 
of ( p.7|) implies 



M, 



^2 u 

^^1 



with 



R 



1 + Egtan/3 + E^ tan/? |hu|^ . 



(2.8) 
(2.9) 

(2.10) 



Notice that ( p.9| ) plays the role of a re-defining (renormalization) condition for the d-quark 
Yukawa couplings, in the process of resumming higher-order radiative corrections. Observe 
also that the matrix R cannot be zero, as this would result in massless d quarks. 

With the help of ( |2.6| ) and ( |2.9| ), we can now express our original Yukawa Lagrangian 
( p.l|) in terms of the mass eigenstates di^R and Ml,_r in a resummed form: 



-C 



Y 



a/2 - 

— dR MdV+R-i 

Vl 



0* 



$r ( Eg 



Eu|h,| 



Vdi + ^^urIi^ul + H.c. 



— ftan/3<l>;* - <^l*] dRMdV^-R-^V dL + — <I>S* rf^j Md 
+ ^IurIi^ul + H.C. 



V2 



(2.11) 



In deriving the last equality in (|2.11|) , we have employed the relation: R ^(Eg + Eu|huP) = 
(1 -R-i)/tan/3. 

It is very illuminating to see how the FCNC part of ( p. 11 ) compares with the literature, 
e.g. with that obtained in Ref . |0 . To this end, let us first assume that 1 + Eg tan (3^0 
and decompose R~^ as follows: 

1 



R-^ 



l-f Eg tan p 



Xfc 



(2.12) 
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where xfc is the diagonal matrix 

Eu|hu|2 tan/3 

yfc = ^ ■ (2.13) 

(1 + Eg tan/3) R ^ ^ 

Making use of the above linear decomposition of R^^ and the unitarity of the CKM matrix 
in (|2.11|) , the FCNC part of our resummed effective Lagrangian reads 



-C^CNC = ^(tan/3$?* - $°*)d,^m,^ (yjx^^yy + OfcK, + KkrcKij^.L, 

(2.14) 

where Xfc'^'*'' ^^e the diagonal entries of xfc ^^nd summation over i,j = d,s,b is understood. 
The term proportional to Xfc gives the top quark contribution, which is the result of [^] 
and subsequent articles |]TU|,0,|T^. However, we should remark here that the frequently- 



used top-quark dominance approximation cannot be justified from considerations based 
only on minimal flavour- violation models |^,|T0|,|1^,[T^. In fact, the other terms in ( p.l4| ) 



and especially the one proportional to Xfc ^'^ charm-quark contribution become 
rather important in the limit 1 -|- Eg tan /3 0. In this limit, the singularity in Xfc i^ 
canceled against the singularities of Xfc Xfc ^ result of the unitarity of V. In this 
context, we should note that the limit 1 -|-Egtan/3 — > is not attainable before the theory 
itself reaches a non-perturbative regime. Requiring that all rf-quark Yukawa couplings are 
perturbative, we can estimate the lower bound, |1 + Egtan/3| > 2.5 x 10~^, for tan/3 = 50.0 
Although |1 -|- Eg tan /3 1 must not vanish in perturbation theory, it can be sufficiently close 
to zero, so that the c-quark contribution becomes competitive with the t-quark one. 

So far, we have assumed that the radiatively-induced matrices Eg and Eu in the effective 
Yukawa Lagrangian (|2.1| ) are proportional to the unity matrix. However, this assumption 
of flavour universality is rather specific. It gets generally invalidated by the mixing of the 
squark generations, the soft trilinear Yukawa couplings and renormalization-group (RG) 
running of the soft SUSY-breaking parameters from the unification to the low-energy scale. 
In this respect, the minimal-fiavour- violation hypothesis, although better motivated, should 
also be viewed as a particular way from minimally departing from universality. 

Nevertheless, given that threshold radiative effects on the up sector are negligible, espe- 
cially for large values of tan /?, it is straightforward to derive the general resummed form for 

^To obtain this lower limit, we simply take the trace of the square of (|]^) and demand that 
Tr|hdp < 3, or equivalently Tr[(R-i)ty|MdP^^^R~^] < 3i;?/2. The latter implies that |1 + Egtan/3| > 
{nis tan (3) /{^/Snit) = 5 x 10~^ tan (3. Finally, it is amusing to notice that if GmEg — and SfteEg < 0, 
a perturbative upper bound on tan/3, tan/3 < l/|Eg|, may be derived beyond the tree level. 
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the Higgs-mediated FCNC effective Lagrangian. Starting from general non-diagonal ma- 
trices Eg and Eu in and following steps very analogous to those from ( p.l|) to (|2.11|) , 
we arrive at the same form as in for the resummed FCNC effective Lagrangian, but 

with Eg, Eu and R replaced by 

(2.15) 



E„ 



R 



E, 



1 + tan/3 Eg + Eu Ihul 



(2.16) 



where the ellipses in (|2.16| ) denote additional (generically sub-dominant) threshold ef- 
fects pO[. Notice that the unitary matrix in ( |2.15|) , which is only constrained by 
the relation |hup = W^^hJ,huW^, introduces additional non-trivial model dependence in 
the matrices Eg, Eu and R. In other words, the presence of reflects the fact that the 
3 X 3-dimensional matrices h|jhu, and m| cannot be diagonalized simultaneously, 
without generating FCNC couplings in other interactions in the MSSM Lagrangian, e.g. 
in the W^-ul-ul and W'-ui-di couplings. Moreover, even in minimal flavour- violating 
scenarios, the 3 x 3-dimensional matrix R may generally contain additional radiative ef- 
fects proportional to W^^h^hdW^ induced by RG running of the squark masses from the 
unification to the soft SUSY-breaking scale Msusy- These contributions can be resummed 
individually by taking appropriately the hermitian square of the modified (|2.9|) and solving 
for W^^hJjhdW^. This last step may involve the use of iterative or other numerical methods. 

In the general case of a non-universal squark sector, the resummed FCNC couplings 
of the Higgs bosons to down- type quarks can always be parameterized in terms of a well- 
defined set of parameters at the electroweak scale. In the weak basis, in which = 1, 
the set of input parameters consists of: (i) the soft squark mass matrices m?^ ^, m| and 
soft Yukawa- coupling matrices Au^d', (ii) the u- and (i-quark masses; (iii) the CKM mixing 
matrix V. 

In our last step in deriving the resummed FCNC effective Lagrangian, we express the 
Higgs fields $5 2 in terms of their mass eigenstates -^1,2, 3 and the neutral would-be Goldstone 
boson in the presence of CP violation Following the conventions of we relate 
the weak to mass eigenstates through the linear transformations: 

1 



V2 
1 

7! 



$? = ^ OuHi + qcos/?G° - sin (3 OsiH, 



02iHi + q sin/3 G° + cos/3 03iHi 



(2.17) 



where Oij is a 3-by-3 orthogonal matrix that accounts for CP-violating Higgs-mixing ef- 
fects IQ. If we substitute the weak Higgs fields by virtue of (|2.171) into ( p.ll|) , we 
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obtain the general resummed effective Lagrangian for the diagonal as well as off-diagonal 
Higgs interactions to d quarks, 



3 



i^H.M' = -l^Y.H^^{^<^ ^H.M' Pl + ^Im' Md'Pii ) d' , (2.18) 
where Pl{r) = [1 - (+)75]/2 and 

cos p V / sm p V COS"' p / tan p 



(2.19) 



The 3 X 3-dimensional matrix R in (|2.19| ), which resums all tan /9-enhanced finite radiative 



effects, is given by ( p.l6|) . Equation ( |2.18|) , along with ( |2.19|) , constitutes the major result 
of the present paper, which will be extensively used in our phenomenological discussions 
in Section 3. 

Finally, let us summarize the most important properties of the resummed effective 
Lagrangian ( ^.181 ): 



i) The FCNC interactions in (|2.18| ) are described by tan^ /3-enhanced terms that are 



proportional to On and O^i and to V^R^^V in (|2.19|) . These tan^ /3-enhanced FCNC 
terms properly take into account resummation,^ non-universality in the squark sector 
and CP-violating effects. 

(ii) The resummation matrix R controls the strength of the Higgs-mediated FCNC ef- 
fects. For instance, if R is proportional to unity, then a kind of a GIM-cancellation 
mechanism becomes operative and the Higgs-boson contributions to all FCNC 



observables vanish identically in this case. Furthermore, as well as the top quark, the 
other two lighter up-type quarks can give significant contributions to FCNC tran- 
sition amplitudes, which are naturally included in ( |2.18| ) through the resummation 
matrix R. 



•^In addition to the non-holomorphic contributions we have been considering here, there are in general 
holomorphic radiative effects on the $5 couphng to d quarks which have an analogous matrix structure, 
i.e. Eg + eu|hup. These additional holomorphic terms are generally small, typically of order 10~^, and 
only slightly modify the form of the matrix R to: R = 1 + (1 + eg + eu|hup)~^(Eg + Eu|huP) tan/3. 
Obviously, such a modification does not alter the general form of the effective couplings in ( ^.19 ) and is 



beyond the one-loop order of our resummation. Therefore, these additional small holomorphic terms can 
be safely neglected. 
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(iii) In the CP-invariant limit of the theory, the effective couphngs g^' j^/ are either pure 
real or pure imaginary numbers. Moreover, in the limit V — > 1, the effective La- 
grangian ( p.l8| ) of the diagonal Higgs-couplings to down quarks is in excellent agree- 
ment with the one presented in . 



(iv) If Mh+ ~ 3 ^ Mhi, one can show that On ~ cos P, O21 ~ sin (3 and O31 ~ 0. In 
this case, the i/i-coupling to d quarks becomes SM like and so iJi-mediated FCNC 
effects are getting suppressed. Instead, the FCNC couplings of the heavy _f/'2,3 bosons 
to (i-type quarks retain their tan /3-enhanced strength in the above kinematic region. 

(v) The one-loop resummed effective Lagrangian ( p.l^ ) captures the major bulk of the 



one-loop radiative effects ||2^ for large values of tan /3, e.g. for tan (3 ^ 40, and for a soft 
SUSY-breaking scale Msusy much higher than the electroweak scale . In addition, 
( |2.18|) is only valid in the limit in which the four-momentum of the d quarks and 



Higgs bosons in the external legs is much smaller than Msusy- This last condition is 
automatically satisfied in our computations of low-energy FCNC observables. 

In the next section, we will study in detail the phenomenological consequences of the 
tan /^-enhanced FCNC effects mediated by Higgs bosons on rare processes and CP asym- 
metries related to the K- and i?-meson systems. 



3 Applications to K- and 5-meson systems 



We shall now analyze the impact of our resummed effective Lagrangian ( |2.18|) for Higgs- 
mediated FCNC interactions on representative K- and 5-meson observables. For compre- 
hensive reviews on K- and i?-meson physics, we refer the reader to [E3,|26|,E7 . 



3.1 AMi^, ck and e'/e 

Our starting point is the effective Hamiltonian for the AS* = 2 interactions, 

Hti=' = ^2MwY. ^^(/^) ^^(/^) ' (3-1) 

i 

where Gp = 1.16639 x 10~^ GeV~^ is the Fermi constant, and Cj(yu) are the scale-dependent 
Wilson coefficients associated to the AS" = 2 quark-dependent operators Qi. Note that the 
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CKM matrix elements in (B.lD have been absorbed into the Wilson coefficients. The AS = 2 



operators Qi may be 


summarized as follows: 






Q^^ = 




gr^ = 


is^^PRd){srPRd), 


Q^^ = 




g^^ = 


{sPLd){sPRd) , 


gf^ = 


{sPLd){sPLd) , 


gSRR ^ 


isPRd){sPRd) , 


gf^ = 


{sa^,PLd){sa^''PLd), 


gf ^ = 


isa^,Pnd)isa''''Pnd) , 


with cr^^ = |[7/.,7i.]- 


Here, much of our discussion and notation follows Ref. 28 



(3.2) 



It now proves convenient to decompose both the K^-K^ mass difference AMk and the 
known CP-violating mixing parameter ex into a SM and a SUSY contribution: 



AMk = Mwr - Mks = AM|^ + AM|u^^ 



ex 



SM I ^SUSY 



'K 



+ e 



K 



To a good approximation, one has 



AM 



SM, SUSY 



2^e{K''\H^J^='\K%M,svsY 



SM 



SUSY exp(i7r/4) 



53m (ir°|ff,^^=2|i^0)s^,_susY. 



(3.3) 

(3.4) 
(3.5) 



V2AMk 

The SUSY contribution to the matrix element {K'^\H^^^'^\K'^)svsy niay be written down 



as 



2|,|9| 



/ f>0\ ttAS=2 I T^Q\ 

I ^cS /SUSY 



'-^F ^^2 _ 172 „ 5 dVLL ( ^VLL , ^VRR 

127r2 



pLR ^LR 



+ PLR^LR 



SRR 
2 



, (3.6) 



where mx = 498 MeV, = 160 MeV and the P's are the next-to-leading order (NLO) 
QCD factors that include the relevant hadronic matrix elements |l28|, ^|3ll , ^[29(] . At the 
scale /i = 2 GeV, they are given by p8| 



jLR 



pVLL ^ Q 25 ^ pLR ^ g ^ pLR ^ 3Q g ^ pSLL ^ _g 3 ^ pSLL ^ _^g g ^ 

On obtaining ([3^71) , we have used the numerical values: ri2 = 0.57 and Bk = 0.85 ± 0.15. 
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From studies in the CP-conserving MSSM with minimal flavour violation [p9ip!0| , [l4| , it is 
known that for large values of tan j3 ^ 40, the dominant contribution to {K^\H^^^'^\K^)svsy 
comes from Higgs-mediated two-loop double penguin (DP) diagrams proportional to 6*2^. 
Within the framework of our large-tan /5-resummed FCNC effective Lagrangian p.l8| ) that 
includes CP violation, the Wilson coefficients due to DP graphs are found to be 



^SLL (DP) ^ _ IGn^mj 9Hisd9Hisd 

' ~ V2GfM^ ^ Ml 

^SRR(DP) _ ISTT^m^ ^ 9Hisd9Hisd 



^LR(DP) _ _ S2n'^msmd 9HiSd9Hisd /o q\ 



where the tan^ /3-enhanced couplings 9^^^ may be evaluated from ( p.l9| ). Note that the 



DP Wilson coefficients in ( p.8|) exhibit a tan"^ /^-dependence and, although being two-loop 
suppressed, they become very significant for large values of tan/? > 40. 

In addition to the aforementioned DP contributions due to Higgs-boson exchange 
graphs, there exist relevant one-loop contributions to {K^\H^^'^\K^)s\}sy at large tan/3: 

(i) the t-H^ box contribution to C^^ of the two-Higgs-doublet-model (2HDM) type, and 

(ii) the one-loop chargino-stop box diagram contributing to Gf^^. The first contribution (i) 
becomes significant, up to 10%, only in the kinematic region Mii± ^ rrit. In this case, to a 
good approximation, Gf^^ may be given by p9|, p!0|] 



^LR(2HDM) _ '^nis'md f^^^-r^ ^2 . 2 n /q n^ 

C2 ~ [Vtsytd) tan 13 . (3.9) 



Note that the light-quark masses contained in (|3.8| ) and (|3.9|) are running and are evalu- 
ated at the top-quark mass scale, i.e. ms{mt) ~ 61 MeV, md{mt) ~ 4 MeV. The second 



contribution (ii) becomes non- negligible only for small values of the /^-parameter [|T0|,|29 
i.e. for < 200 GeV. 

In view of the above discussion, the kinematic parameter range of interest to us is: 

MsusY, /i > m^, tan/3 > 40, (3.10) 

including the case Mh± ^ rrit, for which the Higgs-mediated DP effects can dominate 
the K^-K'^ transition amplitude. Thus, taking also into account the sub-dominant 2IIDM 
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contribution (|3.9| ), formula (|3.6| ) simplifies to 
{K'\ H^/=' \K%ysY = 4.6 X 10-1^ GeV ( 



V2 \(B 



X 



160 MeV / V 0.57 / V 0.85 

30.6 ( + C2"^(2HDM)\ _ /^SLL(DP) ^ ^SRR(DP) . _ ^g^^^^ 



Observe that the Wilson coefficients Ci and C2 contribute with opposite signs to 
SUSY. Based on ( |3.11| ), we will give numerical estimates of AMk and 



ex in Section 4. 

We now turn our attention to the computation of the direct CP-violation parameter 
e'/e in the kaon system, induced by CP-violating Higgs-mediated FCNC interactions. In 
the SM, unlike the Z-penguin graphs, the Higgs-penguin contribution to K ^ -kti , which 
is proportional to the operator 

Qh = {sPLd) J2 (^^) ' (3-12) 

q=u,d,s 

has a suppressed Wilson coefficient proportional to uig/Mfj, where the SM Higgs-boson 
mass Mh is subject into the experimental bound: Mh ~ 114 GeV. One might even think 
of the possibility that the operator Qh in ( |3.12|) , which has enhanced Wilson coefficients 
for q = c,b, mixes with the gluonic and electroweak penguins, as well as with the other 



basis operators in ( |3.2| ). However, as was already pointed out in |Q, this is not the case, 
and so the SM-Higgs penguin effects remain negligible. 

The situation changes drastically in the MSSM with explicit CP violation, since the 
Higgs-boson FCNC couplings to down-type quarks are substantially enhanced by tan^ /3- 
dependent terms, for large values of tan (3. Furthermore, besides the CKM phase, the 
presence of complex soft SUSY-breaking masses with large CP-violating phases may further 
increase the Higgs-boson FCNC effects on e'/e. In fact, we note that soft CP-odd phases 
could even be the only source to account for direct CP violation. 

To reliably estimate the new SUSY effect on e'/ e due to Higgs-boson FCNC interactions, 
we normalize each individual contribution with respect to the dominant SM contribution 
arising from the operator Qq = Yliq=uds^^^R'i)^'i-^Ld) with Wilson coefficient uq, 



VIZ. 



e' /e' 



e (7)6 (^^^ ~'~ ^^^^"^ ~'~ ^sifsY ) ■ (3.13) 
In the MSSM with minimal flavour violation, the non-Higgs SUSY contribution fisusY 
is small |^|. Sizable contributions may be obtained if one relaxes the assumptions of 



universality and CP conservation in the squark sector |37, 38 . 
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Here, we compute a novel contribution to e'/e, namely the quantity ^^su^|y in (|3.13|) , 
which entirely originates from Higgs-boson exchange diagrams in the CP-violating MSSM. 
Based on our resummed FCNC effective Lagrangian ( |2.18| ), we obtain in the zero strong- 
phase approximation 

3 

2=1 q=u,d,! 



oHiggs 
^ ^SUSY 



rris nig 
Ml 



+ 



sd 



ye{7r+7r-\Qe\KO)o 



{n+n-\{sPLd){qj,q)\K°), 
ye{n+n-\Qe\K^)o 



1 {n+n-\{sPLd){qq)\K^)2 
\uj\ ye{7T+7r-\Qe\K^)o 

1 {7i+7i-\{sPLd){q-f,q)\K')2 



\UJ\ 



i=l q=u,d,! 



rrid rriq 
Ml 



(3.14) 



^m{9lM)9kqq [ {7T+n-\{sPRd)iqq)\K^)o 



+ 



A^X^T] 

Hisd) 9H,qq 



ye{7r+7r'\Q,\K^)o 
{n+7r-\{sPRd){q-f,q)\K^)o 



1 {n+7,-\{sPRd){qq)\K')2 
\uj\ ye{n+Tr-\Q,\K<^)o 

1 {n+n-\{sPRd){qj,q)\K% 



\uj\ 



y6{n+n-\Qe\K<^), 



A^X^T] V yG{7r+Tr-\Qe\KO)o 
where the subscripts 0, 2 adhered to the hadronic matrix elements denote the total isospin 
I of the final states, and A'^X^rj is a CKM combination in the Wolfenstein parameterization, 
which has the value A'^X ^r/ ^ 1.3 X 10"^ in the SM. Furthermore, for Aqcd = 325 MeV 
and nis = 150 MeV, the SM Wilson coefficient y^ and the matrix element of take 
on the values ||3^: ye = —0.089 and (7r+7r^|Q6|-^°)o = —0.35 GeV^, respectively. Also, 
experimental analyses suggest the value \uj\ = 0.045, approximately yielding the SM con- 
tribution QsM ~ 1 to e'/e in ( |3.13| ). Finally, the parameters f?^^^ and g^^^ that occur 
in ( |3.14| ) are the diagonal scalar and pseudoscalar couplings of the Hi bosons to u- and 
d-type quarks |jl9[ , whose strengths are normalized to the SM Higgs-boson coupling. These 
reduced ifj-couplings are given by 



9H^dd 
9HiUu 



p 

9Hidd 



2 (9H,dd 



9H,dd ) 



2 ( 9Hidd + 9H^dd ) ' 

02i/sin/3, = -Ogicot/?, 

where we have neglected the small radiative threshold effects in the up sector. 



(3.15) 
(3.16) 



On the experimental side, the latest world average result for 9fJe(e'/e) is [3S 



(1.66 ±0.16) X 10 



-3 



(3.17) 
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at the 1-cr confidence level (CL). In the light of the experimental result ( p. 171) and the 
discussion given above, we may conservatively require that 



< 1 . (3.18) 



The biggest contribution in the sum over quarks in ( p.l4|) comes from the (i-quark and 
exhibits the qualitative scaling behaviour 

^Sy - ^^x^W- (3-19) 



For instance, for tan/3 = 50 and Mh = 200 GeV, (^J^) gives l^g^fy - 0-1 x ^'(1). 



Obviously, such a contribution is, in principle, non- negligible, but very sensitively depends 
on the actual values of the new hadronic matrix elements: 

iQkp)^ = {^^^-\isPL,Rd)iqil,^,)q)\K')j, (3.20) 

with J = 0, 2. A detailed calculation of the hadronic matrix elements {Qs'p)i in ( p.20| ) will 
be given elsewhere. In Section 4, however, we will present numerical estimates of AMx 
and ex within the context of generic soft SUSY-breaking models. 

3.2 AMsg, Bq £^£^ and associated CP asymmetries 

We start our discussion of a set of observables in the S-meson system by first analyzing 
the -Bq--Bq mass difference, AMsg with q = s,d, in the CP-violating MSSM at large tan jS. 
In the applicable limit of equal 5-meson lifetimes, AMb^ may be written as the modulus 
of a sum of a SM and a SUSY term: 

AMsg = 2KBJ|i7f/=2|i?J)sM + (5j|if,^/=^|i?J)susY|, (3.21) 

where the effective AB = 2 Hamiltonian H^^"^ may be obtained from the AS* = 2 
one stated in (|3.1| ), after making the obvious replacements: s h and d ^ q, with 
q = d,s. Proceeding as in Section 3.1, we arrive at analogous closed expressions for the 
SUSY contributions to the AB = 2 transition amplitudes: 
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{B'jH^s^='\B%^SY = 2310 ps 



-1 



265 MeV 



r}B_ 
0.55 



0.88(C2'^(°^) + C2'i^(2HDM)j _ o.52(cf^^(°P) + cfi^(DP))] . (3.22) 

In deriving ( p.22|) , we have also substituted the values determined in [|28|,^^,^,^ for 
the NLO-QCD factors, along with their hadronic matrix elements at the scale /i = 4.2 GeV: 



pLR ^ -0.58, P^^ = 0.88, Pf^^ = -0.52, P^^^ = -1.1. (3.23) 



)LR 



jSLL 



)SLL 



Moreover, the corresponding Wilson coefficients appearing in (|3.22|) may be recovered from 
those in ( p.8|) and ( p.9|) , after performing the quark replacements mentioned above. 

Another observable, which is enhanced at large tan/3, is the pure leptonic decay of B 
mesons [|,|,|,|,|lO|,|n],[T|,|l|,|l4|,||], 5° ^ i+i', with i = fi, t. Neglecting contributions 
proportional to the lighter quark masses the relevant effective Hamiltonian for AB = 1 
FCNC transitions, such as 6 — >• qi'^i~ with q = d,s, is given by 



i7,^«=^ = -2V2GFVtbV,l(^CsOs + CpOp + Cio Oio) , 



where 



(3.24) 



Os 
Op 



O 



10 



167r2 
167r2 
167r2 



mfe (gP^fe) (££) 



m, {qPph) ihd) 



(3.25) 



Employing our resummed FCNC effective Lagrangian (|2.18|) , it is not difficult to compute 
the Wilson coefficients Cs and Cp in the region of large values of tan/3:0 



C, 



3 



2=1 



Ml 



(3.26) 



''Our approach to Higgs-mediated FCNC effects presented here may be extended to consistently account 

for charged-lepton flavour violation in B- meson decays, such as Bs,d ~* £^i'^ |^^, where the effective off- 

s P 

diagonal Higgs-lepton-lepton couplings g^' can be derived by following a methodology very analogous 
to the one described in Section 2. 
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while CiQ = —4.221 is the leading SM contribution. In analogy to (|3.15|) , the reduced scalar 
and pseudoscalar Higgs couplings to charged leptons (7^^^ in (|3.26|) are given by 



9L 



COS P 



p 

9Hd 



tan P O. 



3i 



(3.27) 



where non-holomorphic vertex effects on the leptonic sector have been omitted as being 
negligibly small. 

With the approximations mentioned above, the branching ratio for the meson decay 
to acquires the simple form p[ 



(3.28) 



'tq\ 



'1- 



+ + 2m,Fjp 



where tb^ is the total lifetime of the Bq meson and 



Ba 



■a 



S,P 



r, , ~^,-ri 'A o-^-BgC'lO- (3.29) 

2 ' m;, + TUq 2 

In our numerical estimates, we ignore the contribution from Cio, as being subdominant in 
the region of large tan/3, i.e. for tan (3 40, where all Higgs-particle masses are well below 
the TeV scale. The SM predictions as well as the current experimental bounds pertinent 
to B{B^ —>■ £~^i~) can be read off from Table 1 in 

In the CP-violating MSSM, an equally important class of observables related to 
B{B^^ i^fr) [|l^] is the one probing possible CP asymmetries that can take place in the 
same leptonic 5-meson decays. The leptonic CP asymmetries may shed even light on the 
CP nature of possible new-physics effects, as the SM prediction for these observables turns 
out to be dismally small of order 10"^ This SM result is a consequence of the fact that 
the CP- violating phase in B^-B^-mixing parameter q/p is opposite to the one entering the 
ratio of the amplitudes Al{r'){B^ 

net CP- violating effect on the observable parameter Xl(r) 
cancels out. 

There are two possible time- dependent CP asymmetries associated with leptonic B- 
meson decays that are physically allowed: 

jrdtT{B^,{t)^ltll) - ir dtT{B^,{t) ^ Ipn) 



{q/p){AL{R)/AL{R)) almost 



and^J.p'' , withL 



/~rftr(i?o(t) 



^r^r) 



(3.30) 



A 



R. Under the assumption that q/p is a. pure phase, one finds 



2Xq '^mXq 



CP 



(2 + Xl) + xl\Xq 



A 



2xq '^mXg 



CP 



(2 + x2)|Ag|2 + x: 



(3.31) 
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where Xq = AMb^/Fb, and 



Mt; (VtbVtl\ PiCs + Cp + 2miC,o/{m,MB^j 



In (|3]3|), /?; = (!- 4m,VM|Ji/^ Mf^ is the dispersive part of the matrix element, 

and Cs,p are Wilson coefficients given in (|3.26| ). The maximal value that the leptonic CP 
asymmetries in ( p.3lD can reach is A^p^ = 1 / ^2 + x'^ and is obtained for '^mXg = |Agp. 
From current experimental data |^3|], one may extract the values Xd = 0.76 and Xg ^ 19 at 
the 95% CL, which leads to ^^^(5,) ^ 5% and ^^^(5^) ^ 62%. 

Within the framework of our resummed FCNC effective Lagrangian, we also improve 
earlier calculations of the CP asymmetries by including B^^-B^ mixing effects through 
/\Mf2\ in ( |3.32| ). According to our standard approach of splitting the amplitude into 
a SM and a MSSM part, we obtain for the SM part 



Mt; = Ms^VbBsFI {V,,V;,f Su , (3.33) 



where Su ~ 2.38 is the value of the dominant mf-dependent loop function for a top-pole 
mass mt = 175 GeV. The SUSY contribution to M^*^^^^^'^ = {B^\ H^/='^ \BfisvsY may be 
obtained from ( |3.22| ). 

In the next section, we will present numerical estimates for the K- and i?-meson FCNC 
observables, based on the analytic expressions derived above. 



4 Numerical estimates 

In this section, we shall numerically analyze the impact of the tan^ /^-enhanced FCNC 
interactions on a number of K- and 5-meson observables which were discussed in detail 
in Sections 3.1 and 3.2, such as AM^, e^, AM^^, AM^^, B^ r^r^, Bg jj,'^ ji~ and 
their associated leptonic CP asymmetries. For our illustrations, we consider two generic 
low-energy soft SUSY-breaking scenarios, (A) and (B). 

In scenario (A), the squark masses are taken to be universal and Eg and Eu are pro- 
portional to the unity matrix at the soft SUSY-breaking scale Msusy- The CP-conserving 
version of this scenario has frequently been discussed in the literature within the context 
of minimal flavour- violation models, see e.g. ||^. 

In scenario (B) we assume the existence of a mass hierarchy between the flrst two gener- 
ations of squarks and the third generation, namely the flrst two generations are degenerate 
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and can be much heavier than the third one. In addition, ahhough not mandatory, we 
assume for simphcity that the model-dependent unitary matrix in ( p.l5| ) is such that 
Eg and Eu become diagonal matrices in this scenario. Clearly, in the limit in which all 
squarks are degenerate, scenario (B) coincides with (A). 

In Fig. § we give a schematic representation of the generic mass spectrum that will 
be assumed in our numerical analysis. More explicitly, we fix the charged Higgs-boson 
Mh+ to the value 200 GeV. Since the effect of the gaugino-Higgsino mixing []T3|,|TB| on 
the resummation matrix R can be significantly reduced for m^^ <^ fi, we ignore this 
contribution by considering the relatively low value m^r ~ 2mt <^ Msvsy in our numerical 
estimates, with rrit = 175 GeV. As can be seen from Fig. |^, the third-generation soft squark 
mass m^, the /i-parameter, the gluino mass rrig and the trilinear soft Yukawa coupling Au 
are set, for simplicity reasons, to the common soft SUSY-breaking scale MgusY, which is 
typically taken to be 1 TeV. 

The soft squark masses of the other two generations are assumed to be equal to rriq in 
our generic framework. To account for a possible hierarchical difference between the mass 
scales m^- and m^, we introduce the so-called hierarchy factor p, such that rriq = pmi = 
pMgusY- As has been mentioned above, models of minimal flavour violation correspond to 
scenario (A) with p = 1. As we will see in detail in Sections 4.1 and 4.2, the predictions 
for the K- and 5-meson FCNC observables crucially depend on the values of the hierarchy 
factor p. Equally important modiflcations in the predictions are obtained for different 
values of the soft CP- violating phases (pg = arg {nig) and (pAu = ^^S i^u)- In addition, the 
K- and 5-meson FCNC observables exhibit a non-trivial dependence on the CKM phase 
^CKM, which is varied independently in our figures. 

Although we primarily use tan/? = 50 and Mfj+ = 0.2 TeV as inputs in our numerical 
analysis, approximate predictions for other values of the input parameters may be obtained 
by rescaling the numerical estimates by a factor 

/ tan/3 \" / 0.2 TeV ,^ ^, 

where the integers n and k depend on the FCNC observable O under study. Such a rescaling 
proves to be fairly accurate for tan/5 > 40 and Mh+ ^ 150 GeV, which is the kinematic 
region of our interest. 
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4.1 AMk and {ckI 

The SM effects on AM^ and |ex| were extensively discussed in the hterature [Q, so we 
will not dwell upon this issue here as well. Instead, we assume that the SM explains well 
the experimental results for the above two observables ||43|| : 



AM^'^P = (3.490 ± 0.006) x lO'^^ MeV , (4.2) 
: (2.282 ± 0.017) X 10"^ . (4.3) 



expi 
I 



Given the significant uncertainties in the calculation of hadronic matrix elements, however, 
our approach will be to constrain the soft SUSY-breaking parameters by conservatively 
requiring that AM^^^ and jel^^^^l do not exceed in size the SM predictions.^ 

To start with, we display in Fig. ^ numerical values for the Higgs-boson DP effects on 
AMk and {exl as functions of the gluino phase argm^, where the hierarchy factor p and 
the phase (pAu = arg (At^b) of the soft SUSY-breaking trilinear Yukawa couplings assume 
the discrete values: (p, (pAu) = (1, 0°), (10, 0°), (1, 90°), (10, 90°), (1, 180°), (10, 180°). 
According to our CP-phase conventions |jT9|, fi is always taken to be positive, while 
the CKM phase ^ckm is chosen to its maximal value 90°. The subdominant one-loop 
2HDM contribution coming from W^-H^ box graphs [cf. ( p.9| )] has also been indicated 
by an arrow in Fig. ^. Predictions for Mh+ and tan/5 values other than those shown in 
Fig. may be approximately obtained by multiplying the numerical estimates by a fac- 
tor xo = (tan/3/50)^ x (0.2 TeV/MH+)^. We observe in Fig. || that the resulting values 
for AMl-USY 

can exceed the experimental error in ( |4.2| ) by one order of magnitude, for 
p = 10 and I^Ac/l, \4>g\ ~ 90°. For the same inputs, |ef^^^| takes on values comparable 
to the experimentally measured one ( [4.3| ). Here, we should stress the fact that universal 
squark-mass scenarios corresponding to p = 1 can still predict sizeable effects on e^- This 
non-zero result should be contrasted with the one of the gluino-squark box contributions 



to AMk and |e/^| which do vanish in the limit of strictly degenerate squarks due to a 
SUSY-GIM-cancellation mechanism. 

In our case of Higgs-mediated FCNC observables, however, the situation is slightly 
different. As we have discussed in Section 2, the size of the FCNC effects is encoded 
in the flavour structure of the 3-by-3 resummation matrix R. Since R is diagonal for the 



^Both AMk and jexj place important constraints on the p-?/ plane of the unitarity triangle. The so- 
derived limits can be used to constrain new physics. In this case, a global fit of all the relevant FCNC 
observables to the unitarity triangle might be more appropriate. We intend to address this issue in a 
future work. 
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scenarios (A) and (B) under consideration, we can expand the tan /5-enhanced FCNC terms 
(Vt R-i V)dd' in (|2l8| ) as follows: 



(4.4) 



where d and d' collectively denote all down-type quarks, with d 7^ d' . For the parameters 
adopted in Fig. 0, the quantities R~\ ^ can be simplified further toPl 



R,. 



R^ 



r; 



1 + 



( 37rp2 



+ 



+ 



327rV 
327r2p2 



tan /? 
tan/3 



+ 



327r2 



-i<i>A, 



tan /3 



(4.5) 



Then, from ( [4.5|) , it is easy to see that the off-diagonal elements of V^R^^V increase 
if (pAu^ 4>g = ^-nd so the effective couplings g^^^,, thereby giving rise to enhanced 
predictions. This is a very generic feature which is reflected in Fig. ^ and, as we will see in 
Section 4.2, also holds true for our numerical estimates of 5-meson FCNC observables. 

Neglecting the small Yukawa couplings of the flrst two generations and making use of 
the unitarity of V, we flnd for d ^ d' 



(VtR-^V),,, oc V:dVu' 







\h |2 


_37r ^ 




327r2 



tan/? . 



(4.6) 



If p = 1, the dominant FCNC effect originates from the second term in the square brackets 
of (|]6D, provided |1 + ^ e"*'^^ tan/?| > 5 x 10~^ tan/? (see also footnote 2). If p > 1, then 
gluino corrections become dominant; they are larger by a factor — 3.6. However, 

between the low and high p-regime, there is an intermediate value of p, where (V"''R~^V)drf/ 
does exactly vanish for d 7^ d' , and so the effective couplings 9^^^,- In this case, one has 
R~^ = R~^ = Rj^ in ([4.5|) , implying that R is proportional to the unity matrix. Then, it 
is (V^R~^V)rf(i/ = 0, as a result of a GIM-cancellation mechanism due to the unitarity of 
the CKM matrix. We call such a point in the parameter space GIM operative point. The 
p value, for which the GIM-cancellation mechanism becomes fully operative, may easily be 
determined from (14.61), i.e. 



2 

Pgim 



1 + 



3^ 

327ra;o 



A<t>g-<t>A,j) 



-1 



(4.7) 



For p > 10, the first two equations in (4.5) may be better approximated by replacing p p/\/2- 
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For the MSSM parameter space under study, there is always a GIM-operative value for 
the hierarchy factor p, iff (pAu — (pg = or ivr. For — (pg = 0, we have pciM < 1, 
whereas it is pgim > 1, for cp^u — (pg = iact, the second case is realized in Fig. ^ 

for Pgim ~ 1.22, where AM^^^ and je^^"^! vanish independently of the value of the 
CKM phase. Here, we should emphasize the fact that the value of Pqim does not depend 
on the FCNC observable under consideration and is in excellent agreement with the one 
determined by ( |4.7] ). Because of the above flavour- universal property of Pqim, one may 
even face the very unusual possibility of discovering SUSY at high-energy colliders, without 
accompanying such a discovery with any new-physics signal in low-energy K- and 5-meson 
FCNC observables. 

It is now instructive to gauge the relative size of the different DP-induced Wilson 
coefficients in ( |3.^ ). For simplicity, let us take p = 1. Then, each individual DP-induced 
Wilson coefficient in (|3.8|) may be approximately given by 



Efn^ ^2 _ fit) ^2 /T^*T^ n2 'sp Oli - Q| + 2 ^ OliO si 



Mi 

I Hi 



i=l i=l 



1 



i=l i=l Hi 

where Xfc t-quark dependent entry of the diagonal matrix xfc defined in ( |2.13| ). 

For Mh+ > 180 GeV, CP-violation and Higgs-mixing effects start to decouple from the 
lightest i/i-sector [|I9|. Moreover, in the region tan/3 > 40, the <l>2-component in the H2- 
and ifa-boson mass-eigenstates is suppressed. As a consequence of the latter, we obtain 

where M^^^ = | {Mfj^+M]j^) and the orthogonality of the O matrix has been used. Since it 
is Oil, O31 <^ 1 in the kinematic region of our interest, then on account of (|4.9| ) and ( 4.10|) 



and for maximal CKM phase 5ckm = 90°, the dominant contribution to AM|-^^^ and 
g^SY QQjj^gg from the last DP expression in ( [4.8|) , namely from the Wilson coefficient 
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6*2 in (pl8|), despite the additional suppression factor uid/ms ~ 1/10 with respect to 
^SLL(DP)_ ^^^^ ^ AMfUSY still receives its largest contribution from C2^^^^\ while 
kl^^^l is dominated by the first DP expression in ( |4.8| ), i.e. from the second DP 

expression in ( ^.8| ) is very suppressed with respect to the first one by two powers of the 
ratio md/nis. From Fig. (|^), we also see that in addition to the CKM phase Sqkm, the soft 
SUSY-breaking CP-phases, such as arg (m^) and arg (Au), may also give rise by themselves 
to enhancements of |AM|.^^"*"| and |e^^^| even up to one order of magnitude. Analogous 
remarks and observations also hold true for the 5-meson FCNC observables which are to 
be discussed in the next section. 

Finally, we should comment on the fact that the 2HDM contribution by itself due 
can only give rise to the undetectably small numerical values, 
|AM|P°^| = 5 X 10"^^ MeV and |e|^™| = 5 x 10^^ (indicated by an arrow in the Fig. D, 
for 5cKM = 90°. 



4.2 AMsg, Bq £^£ and associated leptonic CP asymmetries 

In this section, we will present numerical estimates for a number of 5-meson FCNC ob- 
servables, such as the mass difference AM^^, the branching ratio for Bg —>■ and the 
CP asymmetries associated with the 5-meson leptonic decays. The current experimental 
status of these observables is as follows ||43|| : 



and 46 



AMb^ = 0.489 ± 0.008 ps"\ (4.11) 
AMb, > 13.1 ps^S (4.12) 
B{Bs fi^fx-) < 2.0 X 10"^ (4.13) 

B{Bd^T+T-) < 0.015. (4.14) 

Future experiments at an upgraded phase of the Tevatron collider may reach higher sensi- 
tivity to B{B, up to the lO'Mevel g3,|g,|n|]. 

Let us start our discussion by numerically analyzing the 5-meson mass differences 
AM^^ and AMb^- As in the case of the i^-meson observables, we use the same input 
values as those shown in Fig. |^, i.e. Msusy = 1 TeV, Mh+ = 0.2 TeV and tan/3 = 
50. Then, Fig. |^ displays the combined, SM and Higgs-DP, contributions to AMb^ and 
AMb^ as functions of the gluino phase arg(mg), for 6ckm = 90° and different choices 
of hierarchy factor p and (pAu- Note that the SM contributions alone for 5ckm = 90° 
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are displayed by horizontal dashed lines. Even though the SM predictions for ^Mb^ ^ niay 
adequately describe by themselves the experimental values in ( |4.11| ) and (^]12|), they cannot 
yet decisively exclude possible new-physics contributions due to the inherent uncertainties 
in the calculation of hadronic matrix elements such as those induced by SUSY Higgs- 
mediated FCNC interactions. In particular, we observe in Fig. |^ that SM and Higgs-DP 
effects may add constructively or destructively to the mass differences AM^^ ^ . Similar 
features are found in Fig. ^, where the SM and SUSY Higgs-DP contributions to AM^^ 
and AMb^ are plotted versus the hierarchy factor p, for discrete values of 5ckm and (pAu- 
In our SM CKM-phase convention ||4^, unlike the CKM matrix element Vts-, the matrix 
element Vtd is very sensitive to 5ckm values. As a result, the SM predictions for AMb^ 
strongly depend on the selected value of ^ckm, as can be seen from Fig. |. 

In Fig. ^ we exhibit numerical values for the branching ratios B{B^ f^~^f^~) and 
r+r^) as functions for the gluino phase arg(mg), for Msusy = 1 TeV, Mh+ = 
0.2 TeV, tan/3 = 50, and 5ckm = 90°, where p and (j)Au are varied discretely. Since the 
branching ratios are driven by Higgs-penguin effects in the region of large tan /?, predictions 
for other inputs of tan (3 and Mb+ may be easily estimated by rescaling the numerical values 
by a factor xo = (tan/3/50)^ x (0.2 TeV/MH+)^, for tan/3 > 40. Thus, confronting the 
predictions for B{13^ fJ'^p ) with experiment data in ( [4.13|) , combined bounds on the 
tan j3-MH+ plane may be obtained for a given set of soft SUSY-breaking parameters. As we 
see in Fig. ^, these combined bounds become even more restrictive for large gluino phases, 
aigrrig > 90°, in agreement with our discussions in Section 4.1. 

However, there is an additional factor that may crucially affect our predictions for 
the branching ratios of the decays p^p^ and B'^ t^t', namely the hierarchy 

parameter p. As we show in Fig. H, even for the extreme choice of a gluino phase arg rrig = 
180°, B{B^ — > p^ p~) and -23(5° 1^^ 1^^) can get very suppressed for a specific value of 
p in certain soft SUSY-breaking scenarios that can realize a GIM-operative point in their 
parameter space. As we detailed in Section 4.1, this phenomenon occurs for the universal 
value of p = pgim = 1.22, when (pAu ~ 'Pg = ±180°. As can be seen in Fig. |^, the predicted 
values for i?^ /^^/^~ and B^ t^t~ , where (pAu = 0°, confirm the above observation. 

As we have already mentioned, the observables AMb^ and B{B^ — p^p~) exhibit a 
different scaling behaviour with respect to Mh+ and tan (3, through the scaling factor xq 
in ( [4.11 ). Once the above two kinematic parameters are fixed to some input values, the two 
observables AMb^ and B{Bg P^p ) are then rather correlated to each other, since their 
dependences on the /i-parameter and Au are very similar in the minimal flavour-violating 



case p = 1 [|14[- As one can see from Figs. our numerical analysis agrees well with 
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the above result for p = 1. However, we also observe that the correlation is practically lost 
for p > 1, e.g. close to the GIM-operative points, and/or by the inclusion of CP-violating 
effects, as AM^^ and B{Bg p^p") have different dependences on the soft CP-odd phases. 

Let us now investigate the size of the CP asymmetries in the leptonic iJ^-meson decays 
in the CP-violating MSSM; the corresponding CP asymmetries for the Bg meson are ex- 



perimentally constrained to be rather small, less than 5% (see also discussion after p.32| )). 



In Fig. we display numerical values for the CP asymmetries W^p^^^^^^ and A^^p^^^^^^ 
as functions of the gluino phase arg (rrig), for Msusy = 1 TeV, Mh+ = 0.2 TeV, tan/? = 50, 
and 5cKM = 90°. As usual, we independently vary the parameters p and (pAu to take on the 
discrete values p = 1, 10 and (pAu = 0°,90° and 180°. We find that if the B^-B^ mixing 
is consistently taken into account, the typical size of ^^Tp'* '^^'^^^ and ^^Tp'* '^^'^^^ does not 
exceed 0.7% and 3%, respectively. If B^-B^ mixing is not included, the CP asymmetries can 
reach slightly higher values up to 1.2% and 6%, respectively. The apparent reason for the 
smallness of the CP asymmetries is due to the occurrence of an approximate cancellation 
in the sum Cs + Cp in ( |3.32|) at large tan/?, as the muon velocity is /3^ ~ 1. 

Having gained some insight from the above exercise, one may seek alternative ways 
to enhance the di-muon asjTiimetries Aqp ^^'^^'>^^'^'^^\ ^j^jg g^d, the first attempt 
would be to suppress the effect of the B^ — B^ mixing by considering smaller tan /? values, 
e.g. tan/3 < 10. In this intermediate region of tan/5, the above cancellation in the sum 
Cs + Cp does not occur due to non-trivial CP- violating Higgs- mixing effects and so the 
CP asymmetries A^p '^■^'(^)^-^'(^)'* can be significantly increased. To get an idea of the 
magnitude of the CP asymmetries in this case, we consider the so-called CPX scenario 



introduced in to maximize CP- violating effects in the lightest Higgs sector of an effective 
MSSM. In the CPX scenario, the p-parameter and the soft trilinear Yukawa coupling Au 
are set by the relations: p = 4MgusY and Au = 2Msusy- Thus, for Mfj+ = 0.15 TeV, 
Msusy = 1 TeV, tan/5 = 7, p = 10, = 45° and (f)g = 0°, we find that CP-violating 
Higgs-penguin effects can give rise to the CP asymmetries: 

^g]-Mi/^Z) ^ _g%^ ^[fp^^^«^«) ^ -37% (4.15) 

where B{Ba —>■ p^/i^) = 3.6 x 10^^°, which is an order of magnitude larger than the 
SM prediction [^. Another variant of the CPX scenario of equally phenomenological 



importance utilizes the parameters: tan/5 = 10, p = 10, (pAu = 45° and (pg = 90°, with 
^susY = 1 TeV. In this case, we obtain 

^(^^o^/.+mZ) ^ a'cp'^^^^^'^^ ^ 43% (4.16) 
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and B{Bd —>■ yU^yU^) = 4.1 x 10"^, which is two orders of magnitude above the SM pre- 
diction. The above scenario appears to pass all the experimental constraints, including 
those deduced from LEP2 analyses of direct Higgs searches Most interestingly, 

a possible observation of a non-zero CP asymmetry in the leptonic di-muon channel will 
constitute the harbinger for new physics at B factories. At this stage, it is important to 
comment on the fact that the numerical values stated in ( [4.15P and ( [4.16| ) should be viewed 
as crude estimates, since they are obtained entirely on the basis of our resummed FCNC 
effective Lagrangian (|2.18|) at the intermediate tan/? regime. However, in this region of 
tan/3, we expect additional one-loop effects to start getting relevant, such as supersymmet- 
ric Z-penguin and box diagrams. Even though our initial estimates given above appear 
to yield rather encouraging results, a complete study of the leptonic i?-meson branching 
ratios and the respective CP asymmetries for all values of tan (3 would be preferable. 

In the case of r-lepton CP asymmetries, the velocities of the decayed r-leptons (3r 
is roughly 0.5, so one naturally gets an appreciably higher value for the expression 
PtCs + Cp ~ 0.5Cp in ( p.32| ). As a result, larger values for the r-lepton CP asymme- 
tries are expected. Indeed, in Fig. we display numerical predictions for Aqp ^ ^ and 
Aq^ ^ ^ versus the gluino phase arg(mg), for the same values of the input parameters 
as in Fig. ^ Then, the CP asymmetries Aq^ ^ ^ and Aq^ ^ ^ can be as high as 9% 
and 36%, respectively. 

We conclude this section with some general remarks. In addition to the K- and B- 
meson observables we have been studying here, there is a large number of other FCNC 
observables which have to be considered in a combined full-fledged analysis. For example. 



the decay B |l50|,|51| plays a central role in such a global analysis, because it will 

enable us to delineate more accurately the CP-conserving/CP- violating soft SUSY-breaking 
parameter space favoured by low-energy FCNC observables. In this context, constraints on 
CP-violating SUSY phases from the non-observation of electron and neutron electric dipole 
moments (EDMs) should also be implemented. Large CP gluino and stop phases, as the 
ones considered in our analysis, would require either very heavy squarks with masses larger 
than 5-6 TeV or the existence of a cancellation mechanism among the different one-, two- 
and higher-loop EDM contributions Especially, it has been shown recently [^] 

that if the first two generation of squarks are heavier than about 3 TeV, the required degree 
of cancellations does not exceed the 10% level and hence large CP-violating gluino, gaugino 
and third-generation phases are still allowed for wide regions of the MSSM parameter space. 
Finally, in the present numerical analysis, we have concentrated on scenarios that minimally 
depart from the minimal flavour- violation assumption through the presence of diagonal, but 
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non-universal squark masses. In the most general case, however, the squark mass matrices 
and consequently the resummation matrix R of the radiative threshold effects may not be 
diagonal. Such low-energy realizations with off-diagonal soft squark-mass matrices can still 
be treated exactly within the context of our resummed FCNC effective Lagrangian ( |2.18|) , 
by appropriately considering non-trivial quark-squark CKM-like matrices, such as the 3- 
by-3 unitary matrix in ( |2.9| ). 



5 Conclusions 

We have derived the general form for the effective Lagrangian of Higgs-mediated FCNC 
interactions to d-tjpe quarks, where large-tan /3 radiative threshold effects have been re- 
summed consistently (cf. ( |2.18| ) and (|2.19| )). Our resummed FCNC effective Lagrangian is 



free from pathological singularities, which mainly emanate from the top-quark dominance 
hypothesis frequently adopted in the literature, and has been appropriately generalized 
to include effects of non-universality in the squark sector, as well as CP-violation effects 
originating from the CKM-mixing matrix and the complex soft SUSY-breaking masses. 
In particular, our resummed effective Lagrangian can be applied to study Higgs-mediated 
FCNC effects in more general soft SUSY-breaking scenarios, beyond those that have already 
been discussed within the restricted framework of models with minimal flavour violation. 
Also, an approach to resumming radiative threshold effects, very analogous to the one de- 
veloped in Section 2, can straightforwardly be applied to see-saw SUSY models, so as to 
properly describe Higgs-mediated lepton-flavour-violating interactions. 

Within the context of generic soft SUSY-breaking scenarios, we have analyzed a number 
of K- and B-meson observables, such as AMk,b, ^k, e'/^; B{Bs,d and their associ- 



ated leptonic asymmetries which are enhanced by Higgs-boson FCNC interactions for 
large values of tan/3. We have found that the predictions crucially depend on the choice 
of soft CP-violating phases in a given set of soft SUSY-breaking parameters. For example, 
for certain values of the gluino and stop phases, the predictions can reach and even exceed 
the current experimental limits, whereas for other values of the CP-odd phases the FCNC 
effects can be reduced by one or even two orders of magnitude. Most remarkably, we have 
been able to identify configurations in the soft SUSY-breaking parameter space, such as 
Pgim; where a kind of a GIM-cancellation mechanism becomes fully operative [cf. ( [4.7| )] and, 
as a result of the latter, all Higgs-mediated, tan /3-enhanced effects on K- and 5-meson 
FCNC observables are completely absent. 
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Based on our resummed effective Lagrangian, one may now carry over the present anal- 
ysis to a vast numl:)er of otlier K- and i?-meson observables. Evidently, further dedicated 
studies need be performed in this direction. We expect the obtained predictions to affect 
other low- and high-energy observables, such as measurements of electron and neutron 
electric dipole moments and Higgs-boson searches, as well as studies on cosmological elec- 
troweak baryogenesis and dark matter. It would be very interesting to determine to which 
degree the emerging CP-violating MSSM framework with CP-mixed Higgs bosons medi- 
ating tan /5-enhanced interactions to matter could be potentially responsible for all the 
present and future CP-conserving/CP- violating FCNC effects observed in nature. 
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Figure 2: Schematic representation of the SUSY mass spectrum considered in our numerical 
analysis, where rUq and mi denote the masses of the first two and third generations of 
squarks, respectively. The hierarchy factor p, the phase (f)g of the gluino mass, and the 
phases (pAti, of the soft SUSY breaking trilinear couplings, with (t)At — (l>Ab — (/>Au) 
varied independently (see also discussion in the text). 
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Figure 3: SUSY Higgs-DP contributions to ex and AM^ given in units of 10~^2 MeV 
and 10~^, respectively, as functions of the gluino phase arg(m^), for Msusy = 1 TeV, 
Mh+ = 0.2 TeV, tan/5 = 50 and 5ckm = 90°. As is shown above, the different curves 
are obtained for selected values of p and (pAu ■ The size of the 2HDM effect alone on ex is 
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Figure 4: SUSY Higgs-DP contributions to ex and AM^ given in units of 10^^^ MeV and 
10~^, respectively, as functions of the hierarchy factor p, for MgusY = 1 TeV, Mh+ = 
0.2 TeV, tan/5 = 50, and aigrrig = 180°, where the values of 6ckm and (pAu ^^^^ varied 
discretely. 
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gluino phase arg {rrig), for Msusy = 1 TeV, Mh+ = 0.2 TeV, tan/3 = 50 and 5ckm = 90°, 
where the hierarchy factor p and (pAu ^'^c varied independently as shown above. The SM 
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Figure 9: Numerical values for the CP asymmetries A^p ^ ^ and Aqp ^ ^ as func- 
tions of the gluino phase arg (rrig), for Msusy = 1 TeV, Mh+ = 0.2 TeV, tan/3 = 50, and 
^CKM = 90°, where p and (p^u '^'^^ varied discretely. Also shown is the prediction for the 
CP asymmetries without including B^-B^ mixing. 
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Figure 10: Numerical estimates of the CP asymmetries A^p ^ ^ and A^p versus 
the gluino phase axg{mg), for Msusy = 1 TeV, Mh+ = 0.2 TeV, tan/5 = 50, and 5ckm = 
90°, where p and (pAu io,ke discrete values as shown above. 
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